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Abstract 

Beyond their use as numerical tools, quantum trajectories can be ascribed a degree of 
reality in terms of quantum measurement theory. In fact, they arise naturally from considering 
continuous observation of a damped quantum system. A particularly useful form of quantum 
trajectories is as linear (but non-unitary) stochastic Schrodinger equations. In the limit where 
a strong local oscillator is used in the detection, and where the system is not driven, these 
quantum trajectories can be solved. This gives an alternate derivation of the probability 
distributions for completed homodyne and heterodyne detection schemes. It also allows the 
previously intractable problem of real-time adaptive measurements to be treated. The results 
for an analytically soluble example of adaptive phase measurements are presented, and future 
developments discussed. 

PACS: 42.50.Lc, 03.65.Bz, 42.50.Dv 

1 Introduction 

The concept of stochastic Schrodinger equations (sSEs) [1] or quantum trajectories [2] has generated 
an impressive degree of activity in the quantum optics community over the past few years ^ . To 
one outside the field, it may be hard to synthesize all of this activity, not least because there are 
(in my mind) three distinct interpretations of quantum trajectories. I am here classifying authors 
by the ontological status which they ascribe to their quantum trajectories: not real; subjectively 
real; or objectively real. Alternatively, these three interpretations could be stated as quantum 
trajectories being respectively: a mathematical tool; an application of quantum measurement 
theory; or a theory of microscopic state reduction. In the bulk of this paper I am concerned with 
quantum trajectories as a part of quantum measurement theory (that is, with what I am calling the 
subjective interpretation of quantum trajectories). Note that by "measurement theory" I mean the 
body of mathematical techniques used to describe the process of measurement within the standard 
Copenhagen interpretation [3, 4] of quantum mechanics. I am not concerned with metaphysical 
arguments regarding the adequacy of the Copenhagen viewpoint, nor with physical extensions to 
standard quantum theory. However, interpretational issues have been raised in the context of 
quantum trajectories, and the potential confusion regarding the reality of SSEs may obscure the 
results I wish to present. For these reasons, I will begin with a brief review of the interpretations 
of quantum trajectories. I am not aware of a comparison of this sort being in print. 



1.1 Non-Real Quantum Trajectories 

The "non-real" interpretation of quantum trajectories treats the trajectories simply as numerical 
tools for solving the master equation. This is essentially the attitude taken by Dalibard, Castin and 
M0lmer [1] who introduced the concept, although they do mention that their "Monte Carlo Wave- 
function" approach also provides "new physical pictures" . Note the emphasis on wavefunctions in 
their terminology. The value of quantum trajectories for pure states as a computational tool is that 

♦Current (2003) E-mail: H.Wiseman@griffith.edu.au 

^Related ideas have also been developed by workers outside quantum optics. Space limitations prevent me from 
discussing these. 
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it only takes 2N —2 real numbers to store a state vector in an A^— dimensional Hilbert space, com- 
pared to N'^ — 1 real numbers for a state matrix. Of course, this gain is offset by the fact that many 
trajectories arc rcqiiircd to obtain a reliable ensemble average. However, for a large Hilbert space, 
quantum trajectories may offer a real advantage for numerical computation. One very practical 
application is computing the centre-of-mass motion of a fluorescent atom [5, 6, 7, 8, 9]. 

1.2 Subjectively Real Quantum Trajectories 

In the body of this paper, quantmn trajectories arc attributed some degree of reality, in that 
they accurately model the evolution of a system conditioned on continuous observation. This can 
be called "subjective reality" , because the nature of the trajectories depend on the measurement 
scheme chosen by the observer. Once a particular measurement scheme has been chosen, the 
conditioned evolution of the system (as modeled by quantum trajectories) is objectively real in the 
sense that all observers will agree on the results. The subjectivity is present in the sense that it 
is not permissible to make a statement such as "an excited atom in free space will spontaneously 
emit a photon at a random point in time" unless there is a photodetector to record the emitted 
photon. A different detection scheme, such as heterodyne detection, will register a continuous flow 
of radiation from the atom, not a jump. In other words, the stochastic system evolution is only 
real in so far as the system is a subject of observation; it is not real for the system as an object 
in an observerless universe. As well as trajectories arising from different types of measurements, 
trajectories which do not preserve pure states (because of imperfect detection or bath preparation) 
also arise in this interpretation. 

It is the subjective interpretation of quantum trajectories, in addition to their computational 
convenience, which is the motivation for the work by Carmichacl and co-workers [2, 10, 11, 12, 
13, 14]. Most of this work was numerical, but analytic results were obtained for some particular 
systems and particular measurement schemes [12, 13]. Both quantum jump (direct detection) 
and quantum diffusion (homodyne detection) trajectories were used. The measurement theory 
interpretation was also pointed out by Gardiner, Parkins, and ZoUer [15] for the case of direct 
detection (which was formulated much earlier by Srinivas and Davies [16]). 

There are a number of reasons for wishing to consider the measurement interpretation of quan- 
tum trajectories. Firstly, the simulated photocurrents can be used to calculate such quantities as 
the waiting time distribution for photodetections, which are quite difflcult to calculate by tradi- 
tional methods [17]. Another example is the all-or-nothing fluorescence exhibited in the electron 
shelving of a single atom [18, 19]. The technique of linear quantum trajectories recently published 
by Goetsch and Graham [21] has opened the possibility for at least partly analytical solutions of 
measurement theory problems using quantum trajectories. This is the subject of the body of this 
paper, and I will thus defer further discussion until there. 

A second motivation is that subjectively real quantum trajectories give valuable physical in- 
sight into irreversible quantum processes, different from that offered by the master equation. For 
example, a system as simple as a stationary two-level atom shows strikingly different quantum 
trajectory evolution depending on the chosen method of detection, and that this evolution sheds 
light on the respective measurement results [20]. Another example is that coherences which are 
lost in the non-selective (master equation) picture are seen to be merely distributed among dif- 
ferent quantum trajectories [2, 10, 13, 22, 23]. These reconstructed coherences can restore noisy 
oscillatory dynamics to a system which would otherwise exhibit only damped oscillations [24] . 

Thirdly, the measured results can be used to select certain individual systems from the ensemble. 
In this way, it may be possible to see the coherences which are lost in the non-selective ensemble 
evolution [13, 14]. This could be looked upon as a very basic sort of feedback, in which the measured 
result is used to eliminate some members of the ensemble. The use of quantum trajectories in more 
active feedback schemes is a fourth motivation for subjectively real quantum trajectories, which 
has been pursued by myself and Milburn and others [25, 26, 27, 28, 29, 30, 31]. 

1.3 Objectively Real Quantum Trajectories 

The third attitude towards quantum trajectories is a belief in their objective existence, indepen- 
dent of any measurement scheme. This belief is presumably founded on mistrust of the standard 
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interpretations of quantum mechanics, and a desire to replace subjective collapses due to mea- 
surements with objectively real collapses. The roots of this approach are the dynamical state 
reduction models of Gisin [32], Diosi [33], Pearle [34] and Ghirardi- Rimini- Weber [35], which in 
turn hark back to the original objective interpretation of the wavcfunction by Schrodinger [36]. In 
this desire, proponents of objectively real quantum trajectory theory are akin to hidden variable 
theorists, but unlike the latter, they accept the wavefunction as a representation of reality. It is 
only when irreversible processes occur (described by a master equation) that stochasticity enters 
quantum mechanics, giving rise to quantum trajectories. Because the state vector is considered a 
basic; cilcinicint of reality in this interpretation, it only makes sense to consider quantum trajectories 
which preserve pure states. Thus it is that stochastic Schrodinger equations (as opposed to the 
more general stochastic master equations of subjective quantum trajectories) are essential only if 
one regards them as not real or objectively real. 

At least three pairs of authors have published on this theme. The first is Teich and Mahler 
(TM). They published a single paper [37] presenting a theory of quantum jumps in fluorescent 
atoms. It has been shown to be nonphysical by myself and Milburn in reference [20]. The second 
pair, Gisin and Percival (GP), have published many papers [38, 39, 40], also with other authors 
[41], and have received considerable attention. As pointed out by myself and Milburn, the quantum 
trajectory which GP consider can be given an interpretation in terms of unit-efficiency heterodyne 
detection (or ambi-quadrature homodyne detection) of the electromagnetic field radiated by the 
system [20] . However GP attribute an objective reality to their SSE, saying that it "represents the 
evolution of an individual quantum system in interaction with its environment" . By adding ad 
hoc irreversible processes to a model of a three level atom, they are able to induce their equation 
to exhibit jump-like behaviour [41], as seen in electron shelving experiments [18]. In general, 
however, GP cannot associate their trajectories with the detection of photons for the very good 
reason that they actually correspond to heterodyne, not direct detection. A third pair, Breuer 
and Petruccione (BP), have very recently published two papers [42, 43] in which they claim to 
show that the quantum jump unraveling of reference [15] is the SSE for a Markovian open quantum 
system. Although it is in many ways the most natural unraveling, the quantum jump SSE (which 
under some circumstances corresponds to direct detection [15, 20]) does not have any fundamental 
privileged status. BP derive it by choosing a particular basis in which to diagonalize the bath, 
namely its energy eigenstates. 

The basic problem with the objective interpretation of quantum trajectories lies with the origin 
of irreversible evolution. In standard quantum optics, the irreversibility of the master equation for a 
system is an approximation to the exact reversible evolution for the system and its environment [44] . 
None of the above authors try to define in what contexts the master equation evolution is meant to 
be truly irreversible. If there were some fundamental physical mechanism by which irreversibility 
entered the world, for example at some sufficiently large scale, then it is conceivable that one of 
the three schemes outlined above could be correct ^. The quantum measurement "problem" would 
then be solved, as measiirement results would be determined by the irreducibly stochastic quantum 
evolution. Irrespective of this speculation, the models of TM, GP and BP cannot be objectively 
real for the systems to which the authors apply them ^. Irreversible stochastic evolution, in the 
manner of TM, GP or BP, does not take place within an individual atom. The entanglement 
between an atom and its outgoing field is a fact which is not represented by the models of TM, 
GP and BP. This entanglement is most simply shown by the fact that one can choose different 
detection schemes to observe the atom, as discussed in [20]. Gisin and Percival's scheme only gives 
the correct results for heterodyne detection, Breuer and Petruccione's for direct detection, while 
Teich and Mahler's scheme will not give the correct results for any detection scheme [20]. 

1.4 Outline 

Having discussed the non-real and objectively real interpretations of quantum trajectories, the 
remainder of this paper is devoted to quantum trajectories as a part of quantum measurement 
theory. I begin with a discussion of traditional quantum measurement theory, and why this is 
inadequate to describe practical measurements. This leads to an outline of the theory of operations 

^Of course, these are not the only three conceivable dynamical state reduction schemes; two which have gained 
considerable attention are those of GRW [35] and Penrose [45] . 

^These criticism now appears to be accepted by Percival, as evidenced by reference [46]. 
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and effects, and the different ways in which probability can be incorporated in the theory. Sec. 3 
introduces quantum trajectories as an appHcation of quantum measurement theory to continuous 
observation. The description of such trajectories as linear stochastic Schrodinger equations is 
covered in Sec. 4. Sec. 5 presents the general solution for the linear SSE for the case of an undriven 
system with a large local oscillator used in the detection process. This allows a new derivation 
of the operations and effects for completed homodyne and heterodyne measurements. In Sec. 6 I 
apply the results of Sec. 5 to a measurement scheme which has not previously been considered and 
which is probably intractable without linear quantum trajectories. The measurements which I am 
describing here are adaptive measurements, in which the past pliotocurrent is used to modify the 
present conditions of the measurement. Here I discuss the general approach, and present results 
for the simple but non-trivial example of adaptive phase measurements of a field containing at 
most one photon. Sec. 7 concludes. 



2 Quantum Measurement Theory 

Traditional quantum measurement theory is founded on projective measurements. Consider a 
measurement of fixed duration T with a set of N discrete possible measurement results labeled by 
r. With each result r there is an associated projection operator Ilr{T), such that Ilr{T)Ilg{T) = 
Ilr{T)5r.s- If the initial state of the system (assumed pure for simplicity) is then the 

probability for obtaining the result r is 

PriT) = {i;{t)\Ilr{T)m)). (1) 

This can also be written 

Pr{T) = {i;rit + T)\Mt + T)), (2) 
where the tilde denotes an unnormalized state vector 

\Mt + T))=Ur{T)m)). (3) 
Furthermore, the final state of the system is simply this state, normalized if desired 

Although this theory may be fundamentally correct in some sense, it is inadequate to describe 

practical measurements such as those in quantum optics. That is because the experimenter never 
makes a direct measurement on the system of interest. Rather, the system of interest (such as 
an atom) interacts with its environment (the continuum of electromagnetic field modes), and the 
experimenter observes the radiated field. Of course one could argue that the experimenter does 
not observe the radiated field, rather that the field interacts with a photodetector, which triggers a 
current in a circuit, which is coupled to a CRO, which radiates more photons, which interact with 
the experimenter's retina, and so on up the von-Neumann chain [47]. The point is that at some 
stage one has to consider the measurement as being made. If one attempts to apply a projection 
postulate directly to the atom, one will obtain nonsensical predictions. However, it can be shown 
that assuming a projective measurement of the field will yield results negligibly different from 
those assuming a projective measurement at any later stage, because of the rapid decoherence of 
macroscopic objects [48]. For this reason, it is sufficient to consider the field as being measured. 
Because the field has interacted with the system, their quantum states are entangled. If the initial 
state of the field is assumed known (which is not unreasonable in practice, because it is often close 
to the vacuum state), then the projective measurement of the field results in a measurement of the 
atom. Such a measurement however is not proj(X'tive ^. Instead, we need a more general formalism 
to describe such measurements. This is called the theory of operations and effects [50], which I 
will now present in brief. 

Say the initial system state vector is {ipit)) as above, and that the measurement results are again 
labeled by r. Each r is again associated with an operator flr{T), but it need not be a projection 
operator. I will call these measurement operators. The effect [50] for the result r is defined by 

Fr{T) = 0^(T)tO^(T). (5) 
^This is why the original quantum Zeno paradox of Misra and Sudarshan [49] can be avoided in practice. 
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This operator is used to generate the probabiUty for the result via 



Pr{T) = {m\Fr{T)\m)- (6) 

Conservation of probability yields the single restriction which exists for the measurement operators 
[50], namely 

^F,(r) = i. (7) 

r 

The set of all effects {F{r)} constitutes a positive-operator-valued measure (povm) on the space 
of results r [50] . The state of the system given the result r is 

If one were doing only a single measurement, then the conditioned state \ipr) would be irrelevant. 
However one often wishes to consider a sequence of measurements, in which case the conditioned 
system state is vital. In terms of the state matrix p, which allows the possibility of mixed initial 
states, the conditioned state is 

where Pr{T) = Tr[p{t)Fr{T)] and for an arbitrary operator a, 

J[a]p = apa^. (10) 

The superoperator J'[nr{T)] is known as the operation for r [50]. If the measurement were per- 
formed but the result ignored, the final state of the system would be 

p{t + T) = Y,Pr{T)Pr{t + T). (11) 
r 

Consider a hypothetical experiment involving a sequence of measurements for which one knows 

the set of measurement operators {$1^(2^)}. As a theorist, one may be interested in knowing a) 
the probability for obtaining the possible measurement sequences, and b) the state of the system 
during and at the end of the measurement, conditioned on the results obtained. If initially p{t) = 
\'ijj{t)){ijj{t)\ then the above formalism shows that one could obtain this information as follows. 
Start with \^{t)), and generate the probabilities for the results using equation (6). One would 
then choose results at random using these probabilities, and see what the new conditioned state 
of the system is using equation (8). This is the input into the next measurement. The probability 
for each sequence of results to actually occur is simply the probability that it was generated by 
the simulation procedure. The simulation is repeated as may times as necessary to obtain good 
statistics. Call this method A. 

As in the case for projective measurements, this theory can be alternatively stated as follows. 
The final state of the system conditioned on the result r is represented by an unnormalized state 
vector 

\iJ;r{t + T)) = nr{T)m))- (12) 
The probability for this result having been obtained is then 

Pr{T) = {Mt + T)\Mt + T)). (13) 

The spirit of this formulation is that one could simulate the measurement process by randomly 
generating a result r with probability 1/iV, then calculating the appropriate unnormalized condi- 
tioned state. This is the input into the next measurement step. The actual probability assigned to 
each sequence of measurement results is the final norm squared of the state vector (13), where here 
T would be the total time for the sequence of measurements. This simulation would be less efficient 
because all measurement results are given equal computational time, even those with vanishingly 
small probabilities. However it would have the advantage that probabilities need to be calculated 
only once (at the end). Call this method B. 
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In fact, it is possible to formulate the theory in any number of ways which include these two 
as special cases. Define an (in general unnormalized) conditioned state vector by 



Here the Ar{T) are positive numbers chosen such that 

Y.K{T) = 1. (15) 

r 

They can thus be interpreted as probabilities in some sense, and I will call A^. the ostensible 
probability for the result r. The actual probability for getting result r is 

Pr{T) = K{T){Mt + T)\Mt + T)). (16) 

If the Ar(T) are chosen such that Aj.(T) happens to equal {tp{t)\Fr{T)\-ip{t)) , then this formulation 
is equivalent to the method A given above, with \'il}r{t + T)) = \xl}r{t + T)). If they are chosen 
such that Kr{T) = 1/N, then this formulation is equivalent to the method B, with Itprit + T)) = 
{■iprit + T))y/N. In general, one could do the simulation of the measurement by choosing the results 
r with probability Ar{T), then calculating the conditioned state vector |V5r(t + T)). Each result 
would then be weighted by {'ipr{t + T)\'ipr{t + T)), where again T would be the total time. Because 
of the probability of generation, this weighting gives the correct overall weighting to each member 
of the ensemble as given in equation (16). Call this method C. 

This flexibility of formulation has been known for some time (see Goetsch and Graham [21] and 
references therein). As presented here it must appear rather pointless. However we are about to 
consider quantum trajectories, as an application of measurement theory for continuously (in time) 
monitored systems. In this context, it turns out that the ability to separate out the probability 
of each possible result into two components (a probability of generation Ar{T) and a weighting 
(ipr(t + T)|'!/5r(t + T))), will be of enormous use in solving the quantum trajectories. 

Before turning to the topic of quantum trajectories, there is one more useful fact which can 
be explained in terms of the general formalism. Recall that the final state of the system if all 
measurement results are ignored is given by equation (11). This can also be written as 

p{t + T)^Y.'^i^r]p{t). (17) 
r 

Now if we define a unitary rearrangement of of the measurement operators Or- by 

O; = ^ Ur,sfts, (18) 
s 

where [/ is a c— number matrix satisfying Ur,sU*^q = 5s,q. Then the unconditioned final state 
under the new measurement operators {ilj.} is 

^n;p(t)n;^ = Y.Ur,sfisPmlu:^, (19) 

r 'r',s.q 

= Y,s,,gn,p{t)nl = J2nsp{t)nl (20) 

s,q s 

which is the same as that under the old measurement operators {Or} • That is to say, one 
can unitarily rearrange the measurement operators without changing the non-selective system 
evolution. 



3 Quantum trajectories 

One case of the above measurement theory which is of great importance is continuous observation, 
with an infinitesimal measurement interval T = dt. The form of such a measurement theory 
is restricted by the requirement that the non-selective evolution generated by the measurements 
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should give a valid evolution equation for the state of the system. Explicitly, this non-selective 
evolution is 

p{t + dt) = Y^ nr{dt)p{t)Ql.{dt). (21) 

r 

This is obviously Markovian [51], with the increment in the system state from time t to time t + dt 
depending only on the state of the system at time t. A Markovian evolution equation for the state 
matrix of a system is called a master equation. Physically, it can be dcrivc^d from the interaction 
of the system with its environment, if the environment can be treated as a bath. That is to say, it 
is necessary for the environment to irreversibly carry away information from the system. Alterna- 
tively, one can view the environment as a continuous stream of independently prepared apparatuses 
which are wheeled up to the system, interact for an infinitesimal time dt, and are wheeled away 
again to be measured. Although this sounds unrealistic, the electric-dipole interaction between a 
system and the electromagnetic vacuum behaves like this to a good approximation [15, 31]. 
The most general form of master equation is the so-called Lindblad form [52, 44] 

p=-i[H,p]+Y,1^[c^]p- (22) 

Here H is an Hermitian operator and V is a superoperator taking an operator argument (enclosed 
in square brackets). It is defined for an arbitrary operator r by 

V[r]p = rpr^ - i(rVp -|- prV). (23) 

The operators are completely arbitrary. 

Consider now a special case where there is but one Lindblad operator c for the system, so that 

p = -i[H,p]+V[c]p = Cp. (24) 

Then it is possible to write the master equation as the nonselective evolution due to a measurement 
[as in equation (21)] with just two measurement operators Q,r{dt) by choosing 

Qi{dt) = \/dic, (25) 

Qo{dt) = l-{iH+^c^c)dt. (26) 

This is because 

p{t + dt)= J[^r{dt)]p{t) = (1 + Cdt)p{t), (27) 

r=0,l 

where C is as given in equation (24). 

With the measurement operators Q.o{dt)^rti{dt), it is evident that the measurement record 
will be a point process [53], as I will now explain. For almost all infinitesimal time intervals, the 
measurement result is r = 0, because Pa{dt) = 1 — 0{dt). The result r = is thus regarded as 
a null result. In the case of no result, the system state changes infinitesimally, but not unitarily, 
via the operator Jlo. At randomly determined (but not necessarily Poisson distributed) times, 
there is a result r = 1, which I will call a detection. When this occurs, the system undergoes a 
finite evolution induced by the operator fii. This change can validly be called a quantum jump, 
although it must be remembered that it represents a sudden change in the observer's knowledge, 
not an objective physical event as in Bohr's original conception [54]. Real measurements which 
correspond approximately to this ideal measurement model are made routinely in experimental 
quantum optics. If c is the lowering operator for the dipole of the quantum system, multiplied by 
the square root of the damping rate, then this theory describes the system evolution in terms of 
photodetections [31]. 

Now these measurement operators can be unitarily rearranged by the matrix 



Uw Uii ) \ -j*Vdi l-^\j\^dt 



(28) 



This is equivalent to the transformation 

c^c + j; H ^ H -iU-i*c--fc^), (29) 
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under which the master equation (24) is of course invariant. This transformation preserves the 
property that the collapse operator Cli is infinitesimal, and the smooth evolution operator Oq is 
close to unity. Explicitly, 

Oi(dt) = \/rfi(c + 7), (30) 
Qo{dt) = [i/f+i(c7*-ct7) + i(ct+7*)(c + 7)] . (31) 

Physically, it amounts to adding a coherent field to the output radiation field before it is detected 
by the system. This can be done using a low-reflectance beam-splitter. In the limit |7|^ ^ (c^c), 
this procedure is known as homodyne detection. 

The specification of the measurement operators is all that is necessary to define the measure- 
ment. However, it is useful to give a more explicit description of the evolution of the system under 
this scheme. This is done by introducing a stochastic increment dN[t) which equals one if the 
result 1 is obtained in the time interval {t, t + dt) and otherwise. Formally, dN{t) is defined by 

dN{tf = dN{t), (32) 
E[dN{t)] = {n\idt)ni{dt)) =dt{iPit)\ic + j)\c + j)\^it)), (33) 

where a classical expectation value is denoted by E and the quantum expectation value by angle 
brackets. Prom the measurement operators (30,31), the stochastic evolution equation for the 
conditioned state vector is 



dN(t) 



■7 



\m)- (34) 



d\m) 



This gives the quantum trajectory for the system explicitly as a non-linear stochastic Schrodinger 
equation. Putting 7 = gives the form generally used for numerical solutions of the master 
equation. 



4 Linear Quantum Trajectories 

The evolution described by equation (34) corresponds to the method A of Sec. 2. That is to 
say, the results dN{t) are generated with the distribution they would have in actuality. No other 
weighting is necessary. In this case, method B would be outstandingly inefficient. It would imply 
either generating or not generating a photodetection with equal probability for every instant of 
time {E[dN] = i). The norms of the vast majority of states would rapidly become vanishingly 
small, and their contribution to the ensemble would be negligible. The only states which would 
significantly contribute to the ensemble would be those in which photodetections were sparse, as 
they would be in actuality. However, there is a half-way house, corresponding to method C of 
Sec. 2. This would be achieved by keeping E[dN] = 0{dt). Specifically, a convenient choice is 

Ki{dt) =E[dN] = \-i\^dt, (35) 

so that by necessity h.Q{dt) = 1 — |7pdt. Then the unnormalized conditioned system state would 
evolve via 

m^dt)) = '^^W^ = (^+-)m)) (36) 



VAi(dt) V 7 

\Mt + dt)) = ^}2(^}^ = [l-dt{^c^c + cr + iH)]m)} (37) 

The ostensible statistics for N{t) are those of a Poisson process. The actual probability for a 
detection varies in time according to 

Pi(rft) = Ai(dO(^i(i + rfi)IV^i(* + c«t)) =«^i(^(i)l(c + 7)^(c + 7)l^(t)), (38) 
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as before. 

Writing out this evolution explicitly, we have 



dim) 



dN{t)^ - dt [iH + ic^c + 7*c) 



m))- (39) 



This is considerably simpler than equation (34), and, more importantly, it is linear in 
This has been achieved by transferring some of the information about the probability of each 
measurement record into the norm of \'4){t)). The choice of |7p(it for Ai(dt) is obviously not 
a good one if the local oscillator amplitude 7 vanishes. In the remainder of this paper I am 
concerned with the opposite limit, I7I 00. Consider a short time (5t <C 1 in which the system 
changes negligibly, but the number of detections 5N{t) ^ \^\'^5t is very large. This allows the 
Poisson process 5N{t) to be approximated by a Gaussian process: 

5N{t):^\-f\^5t+\-f\5W{t), (40) 

where 5W{f?) is a Gaussian random variable of zero mean and variance St. Since this time is 
infinitesimal as far as the system is concerned, we can substitute this into equation (39) to obtain 

d\i){t)) = [dW{t)e-''^c - dt (iH + ^c^'c)] \%jj{t)), (41) 

where $ = arg7 and dW{t) can be treated as an infinitesimal Wiener increment [51] satisfying 

E[dW{t)] = , dW{tf = dt. (42) 

This equation has been previously derived by Goetsch and Graham, but in a different way. Their 
method involved considering ideal measurements of one quadrature of the radiated field, without 
involving photodetection or a local oscillator. As there is no known physical mechanism for under- 
taking such an ideal measurement, I believe that it is valuable to be able to derive the linear SSE 
(41) as a limit of a physical detection scheme. 



5 Solving the Linear Quantum Trajectory 

It is now convenient to introduce a new notation. Let the recorded photocurrent be defined 

m = lim «^ = (43) 

and let the complete record of the photocurrent from time to just before time t be denoted I[o,t) . 
Then as a reminder that \'ip(t)) is conditioned on the photocurrent, denote it by l-^;). Goetsch 
and Graham [21] have solved the linear SSE (41) for the case of homodyne detection, where the 
local oscillator phase $ is held constant, and for the case of heterodyne detection where it varies 
rapidly: $(i) = $0 + ^t. This was done for some particularly simple internal Hamiltonians H. 
Here I show that a solution may be found in a similar manner for any variation of $ with time. 
I consider only the case H = 0, which amounts to a freely damped cavity. That is because I am 
interested in describing a single measurement of some property of the initial state of the cavity, in 
which all of the light is allowed to leak out and be detected. 

Consider an initial coherent state I^/jq) = |a) where a\a) = a\a), where a = c is the annihilation 
operator for the damped mode, obeying [a, a^] = 1. Now coherent states are unique with respect to 
the damping master equation (24) with H — 0, in that they are the only states which remain pure 
under that evolution equation. Specifically, they have the solution \ipt) = |ae~*/^). This means 
that the solution to the linear SSE (39) must be of the form 

\i,l)=i,l\ae-*/'). (44) 

Now using the Ito calculus, equation (41) can also be written as 

IViJ+rft) =exp(-ia^adi)exp(dV^(t)e-^*(*)a- ie-2^*Wa2(it)|V;i). (45) 
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Now using the fact that 



exp{-^a^at)\a) = exp ( - - e"*))|ae"*), (46) 

one has 

VSJ+dt = exp ( - i|a|2e-*dt - ^e-^"^'^'^ a^e-'dt + er"^^''' ae-'/^I{t)dt)^j\, (47) 
where it is to be remembered that dW{t) = I{t)dt. The total solution satisiying = 1 is thus 

^\ = exp ( - i|a|2(l - e-*) - ^a^ f e'^'^^'^e-'ds + a f e-'*(^)/(s)ds) . (48) 

Jo Jo 

Now it follows from this that the solution for the initial coherent state |V'o) = \ct) can also be 
written 

li^l) = ex^p{-^a^at)exp{^S*it)a^ + R*{t)a)\i;o), (49) 

where ^ ^ 

R{t) = [ e'^^'h-'/^I{s)ds ; S{t) = - f e^i^We-'^ds. (50) 
Jo Jo 

That is to say, the solution at time t depends not on the full photocurrent record T-[o,t) (which is 
a continuous infinity of real numbers), but rather on that record only through the two complex 
integrals (50) ^. Furthermore, because equation (39) is a linear equation by construction, the 
solution (49) is true for arbitrary initial conditions, not just for an initially coherent state. Thus any 
state conditioned on photodetection in the large local oscillator limit depends on the photocurrent 
only by the two integrals R{t) and S{t), regardless of how the local oscillator phase varies. 

The measurement is complete in the limit t — > oo, when there is no light in the cavity and 
hence no information left to be obtained. In that limit, lim(_>oo cxp(— otoi) = I0)(0|, so that the 
final conditioned state |0) is obviously of no interest. What is of interest is the probability for 
having obtained the photocurrent record i[o,t) ■ As we have seen above, the only relevant aspects 
of this record are R{oo), which I will denote A, and 5(00), which I will denote B. Prom Sec. 3, 
the probability for obtaining the results A and B is 

PiA,B)=Po{A,B){^^^\^^''), (51) 

where |^^'^) is the unnormalized conditioned state of equation (49) and Po{A, B) is the ostensible 
distribution for A and B, given by 

Po{A,B) = J dI^o,^^Po{l[o,oc))S^^\A-R{cx))6^^\B-Si^)). (52) 

Here i?(oo) and 5(oo) are the functionals of I[o,oc) defined above, and Po(l[o,oo)) is the ostensible 
distribution for I[o,oo)) which as seen above, is simply that of Gaussian white noise. 
From equation (49) we can rewrite the probability for obtaining results A, B as 

PiA,B)=Tr[F{A,B)po], (53) 

where po is the initial system state (which is allowed to be mixed), and F{A,B) is the effect for 
the results A, B defined by 

FiA, B) = Po{A, B)Ct\A, B)n{A, B), (54) 

where the (unnormalized) measurement operators are 

n{A,B) = \0){0\exp{^B*a'^ +A*a). (55) 

This is a completely general result describing the sorts of measurements which can be made on 
an optical mode by photodetection in the large local oscillator limit. From it can be derived the 
standard results for homodyne and heterodyne detection (as shown below), and also some more 

®It might be thought that S{t) does not even depend on the record I[o,t)i but it might if the local oscillator phase 
is adjusted according to I[o,t) , so as to make an adaptive measurement as will be considered below. 
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interesting results for adaptive measurements (as I will discuss in Sec. 6). At finite times one has 
less than maximal information about the system and as a result the effect cannot be factorized in 
terms of measurement operators. This is a minor generalization of the theory of operations and 
effects as stated in Sec. 2 [44]; the same condition (7) for the effects applies. Explicitly, equation 
(49) shows that the povm for the measurement up to time t is 

Ft{R{t),S{t)) = Po{R{t),S{t))exp{^S{t)a'<^ + R{t)a'')exp{-a^at)exp{^S*{t)a^ + R*{t)a), (56) 
where Po(-R(i)i '^(*)) defined analogously to (52). 



5.1 Homodyne Detection 

For homodyne detection of a given quadrature of the system, one desires the local oscillator phase 
$ to be constant. This means that B is simply given by — e^** Jg°° e~*dt = — e^'*. Thus it can be 
ignored as a measurement result, and we need consider only 

POO 

A = e'^ e-*'^dW{t). (57) 
Jo 

Since dW{t) is ostensibly white noise, it follows from equation (50) that we can write A = e**X 
where the ostensible distribution for X is 

Po{X)dX = exp(-XV2)dX. (58) 
\/27r 

Thus the effect for such a measurement is 

F{X) = exp(-XV2)^(X)|0)(0|^t(x), (59) 

where E{X) is an operator defined by 

E{X) = exp ( - ie''*a^' + Xe'*a+) (60) 

Now using some simple boson operator algebra [55], it is easy to show that 

(ae-** + a)e'^)E{X)\Q) = E{X){ae-"^ + X)\0) = XE{X)\0). (61) 

That is to say, E{X)\0) is an unnormalized eigenstate of of ae~** + a^e** with eigenvalue X. The 
other factor in equation (59) supplies the normalization, so we can write 

F{X) = \X){X\, (62) 

where \X) is an eigenstate of ae~** +a'l'e** satisfying the delta-function normalization {X\X') = 
6{X — X'). This is as has been known for a long time; that homodyne detection can measure a 
quadrature of the field arbitrarily accurately (see reference [56] for a "traditional" derivation using 
quantum Langevin equations). Specific cases of this result have been derived using the theory of 
quantum trajectories (see for example reference [13]). However, as far as I know, this is the first 
time it has been derived as a completely general result using this theory. 



5.2 Heterodyne Detection 

Under heterodyne detection, the local oscillator phase $(f) varies rapidly at a rate A 1. This 
means that B averages over time to zero. The remaining measurement result A is defined by 



poo 

Jo 



(63) 



It is easy to verify that in the limit A » 1, this has the ostensible distribution 



Po{A)d''A = - exp(-]^]2)d2A (64) 
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Thus the effect for a completed heterodyne measurement is 



F{A) = - exp{-\Af) exp(aU)|0)(0| exp{aA*). 



(65) 



Using the techniques of reference [55], it is again simple to show 



exp(-|^|V2 + aU) 



exp(a'^^ - aA*) exp(a^*), 



(66) 



and thus that 



exp(-|A|V2 + aU)|0) = cxp(aU- 
where this last state is the normalized coherent state a\A) 



aA*)\0) = 

= A\A). Thus we have 



(67) 



F{A) = hA){A\. 



(68) 



In this case the factor of 7r~^ is necessary because the coherent states are overcomplete. This 
result, that the distribution function for completed heterodyne measurements is the Q function 
Tr~^ {A\p\A), has also been known for some time (and again there is a succinct derivation in reference 
[56]). It has not previously been derived using quantum trajectories as here. 

6 Adaptive Measurements 

In the two examples considered above, homodyne and heterodyne measurement, the quantum 
trajectory method for finding the finite time effect simply confirmed already known results. A 
type of optical measurement for which the effects have never before been calciilated, and for which 
the quantum trajectory method is perhaps the only tractable approach, is continuously adaptive 
measurements. These are measurements in which the results so far (the photocurrent record I[o,f) 
up to time t) arc used to alter the imitary matrix which determines the measurement operators 
(28) at time t. Note that, as expressed in equation (20), this does not change the overall evolution 
of the system as given by the master equation (24), and so is distinct from feedback onto the system 
dynamics as explored in Refs. [25, 26, 27]. In the infinite local oscillator amplitude limit being 
considered here, there is only one parameter which determines (28); that is the local oscillator 
phase The feedback-control of this phase is just what is required for a measurement which is 
of particular interest: that of the optical phase of the system. 

A standard measurement of the phase of an optical system consists of making a completed 
heterodyne measurement yielding the result A as above, and then discarding the amplitude infor- 
mation in A. That is to say, the effect for such a standard phase measurement is 



where A = [Ale*"^. Because it discards information, this effect cannot be factorized into measure- 
ment operators as in equation (5). In the language of reference [57], it is an unsharp measurement 
of phase. By contrast, an ideal measurement of phase [58, 59, 60] has the effect 



defined in terms of the unnormalized phase eigenstates [</)). Ideal measurements are the optimal 
way to extract information which has been encoded in the phase of the field [61]. 

Standard measurements of the phase are suboptimal because they measure both the phase and 
amplitude quadratures of the field, but use only the former information. It is intuitively obvious 
that one could make a better measurement (i.e. closer to ideal) by concentrating on measuring 
the phase quadrature. However assuming that one knows the phase quadrature before making 
the measurement is not in the spirit of a true phase measurement. Rather, one would have to 
estimate the phase quadrature during the course of the measurement, from the photocurrent so 
far. This is the essence of the adaptive measurement of phase. Precisely what algorithm one used 




(69) 




OO 



(70) 



n=0 
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to estimate the phase quadrature would depend on the knowledge one had about the initial state 
of the system. Nevertheless, the analysis given above shows that whatever algorithm one used, 
the estimated phase at time t, ipi should depend on the photocurrent I[o,t) only through the two 
complex numbers S{t) and R{t). In order to measure the phase quadrature, the local oscillator 
phase should therefore be set equal to 

m-^^t[R{t),S{t)]+'K/2. (71) 

When the measurement is over (at t = oo), one's estimate for the original phase of the system is 

Exactly how (pt should be chosen as a function of R(t) and S{t) is a complicated issue. As stated 
above, it depends on the knowledge the observer has regarding the preparation of the system. This 
is evident in figure 1, which shows the Wigner function of the effect Ft(^R{t), S{t)) (56) for some 
typical values of R and S which might occur in an adaptive phase measurement. The position and 
shape of this positive operator in phase space {q = a + a^ andp = —ia + ia') indicates the inherent 
uncertainties in the results of the single measurement up to time t yielding results R and S. As 
this figure shows, the state Ft(^R{t), 5'(t)) is a Gaussian state. As t — > cx3 it becomes a minimum 
uncertainty squeezed state. For heterodyne detection this reduces to a coherent state of complex 
amplitude A. For homodyne detection it is an eigenstate of p cos ^ + g' sin 0, where 6 is in this case 
equal to the chosen local oscillator phase <I>. In general 6 = arg(5)/2 is a measurement result which 
changes over the course of the adaptive measurement. The displacements x and y are defined in 
terms of the measurement results R and S by 

2jR(i?e-^^) . 2^{Re-^^) 
^- l_e-t-|5| ' ^ - l-e-* + |5r ^ ' 

The variances in the x and y quadrature are 

Evidently, the phase of this state is not unambiguously defined in general. One reasonable estimate 
would be the phase of the centre of the Wigner function, which would be 

<ft = Ot + a.rcta.n{yt/xt) = arg [R{1 - e"*) + SR*] . (74) 

However, depending on one's initial knowledge of the state of the system, a different estimate may 
be more appropriate. 

There is one example for which a different estimate is definitely more appropriate, and for 
which there is no ambiguity about the estimated phase of the system. This is when there is at 
most one photon in the system (and the observer knows this). First, I wish to present the results 
in this case for ideal and standard phase measurements. Projecting the optimal mode effects into 
the subspace spanned by {|0), |1)} gives the ideal phase measurement effect 

fideal(<^) = ^|<^)(<^| , !</-) = |0) + e'^|l), (75) 



and the marginal heterodyne effect 



/TT 



i^std(<^) = ^^ideal(<^) + (^l-^j^- C^^) 

That is to say, the standard technique has an efficiency of \/7r/2 ~ 88%, in the sense that the 
same POVM would arise from a ideal measurement which worked 88% of the time, and which gave 
a random answer on the interval [0, 27r) the other 22% of the time. 

I have considered the adaptive phase measurement of a field with at most one photon in detail 
in reference [62]. The lack of ambiguity in the estimated phase quadrature (pt is due to the fact that 
with at most one photon in the field, the measurement result S{t) is irrelevant to the estimation 
of its state, which is obvious from equation (49). The estimated phase quadrature in this case is 
ift = argi?(t). Substituting this into equation (71) leads to an analytic solution for Po(^): 

Po{A)(fA = 7r-i(5(|A|2 - l)(fA. (77) 
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In other words, A can contain only phase information, not amphtude information. This is precisely 
what is desired for a phase measurement, and indeed the povm for this adaptive measurement is 



Here I have again used the result (54) for an optical mode, projected into the lowest two levels 
{|0), Thus for any field with at most one photon, a sinipk; adaptive measurement scheme can 
produce an optimal measurement of the phase, which is significantly better than the result with 
no feedback (76). 

7 Conclusion 

In this paper I have investigated the reciprocal links between quantum trajectories and quantum 
measurement theory. The first link is that one interpretation of quantum trajectories is as a special 
case of quantum measurement theory, that of continuous observation. Here each measurement has 
a duration which is infinitesimal on the time scale of the non-trivial system evolution, and is 
described in terms of operations and effects. This leads naturally to discrete measurement results 
(photodetections). The equations are actually easy to solve in the limit of a very large local 
oscillator, in which case the discrete photodetections are replaced by a noisy photocurrent. This 
is the limit which was considered in detail here. 

The second link is that by solving the quantum trajectories over a finite time one can obtain 
the operations and eff'ects for measurements of finite duration. These results would be difiicult 
to obtain were it not for a formulation of quantum trajectories as linear stochastic Schrodinger 
equations due to Goetsch and Graham which I have rederived here in a new way. For the case where 
the the system is not driven, but simply allowed to radiatively damp, the finite measurement time 
is naturally extended to infinity in order to detect all of the radiated light. Then it is possible to 
rederive the standard results for homodyne and heterodyne detection, and also to derive new results 
which would be very difficult to derive by other means. These new results are those pertaining to 
adaptive measurements, in which, for example, the photocurrent up to time t is used to control the 
local oscillator phase at time t. Here I presented the results for one case which can be completely 
solved analytically: adaptive phase measurements on a system with an upper photon number 
bound of one. 

The topic of real-time adaptive measurements begins rather than ends with phase measure- 
ments. Feedback-control of measurements can only improve the quality with which properties 
are measurable compared with standard non-adaptive techniques. For measuring quantities other 
than phase, one would not wish to work in the infinite local oscillator amplitude limit. Rather, 
one would have a finite local oscillator, in which case there would be two parameters to be con- 
trolled: its amplitude and phase. This is obviously a more difficult problem, and it may be that 
only a numerical solution using stochastic Schrodinger equations would be possible. Even using 
numerical techniques, it is not obvious what the general algorithm would be for optimizing such 
measurements, nor even how optimization should be defined. These are problems for future work, 
continuing the fruitful interaction between formal quantum measurement theory and its principal 
realm of application in quantum optics. 
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Figure 1: One-standard-deviation contour for the Wigner function of the effect Ft{R,S) (56) for 
typical values of R and S which would occur in the adaptive measurement of the phase of a state 

with a moderate coherent amplitude. Here q = a + and p = —ia + ia^ , x and y are defined by 
equation (72), and the variances in these quadratures by equation (73). 
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